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A DIGITAL PROGRAM FOR CALCULATING THE INTERACTION 
BETWEEN FLEXIBLE STRUCTURES, UNSTEADY AERODYNAMICS AND ACTIVE CONTROLS 

I 

Ell wood L. Peele and William M. Adams, Jr. 


INTRODUCTORY REMARKS 


Active controls technology, (ACT), applied to airplanes involves the 
interaction between a flexible structure, the aerodynamic forces acting 
thereon, and automatically controlled aerodynamic surfaces. The aerodynamic 
surfaces are designed to alter the amplitude, phase and/or spatial distribution 
of these forces. The objectives are to achieve satisfactory rigid body and 
aeroelastic behavior using control surfaces in place of added structural 
material or less efficient planform geometry. There are five major appli- 
cations of ACT ; augmented rigid body stability, augmented flexible body 
stability, reduction of rigid body gust response, reduction of structural 
loads due to gust, and the reduction of wing loads due to maneuvers. All of 
these tasks require the control of one or more aerodynamic surfaces responding 
to a combination of one or more rigid or flexible airframe responses or gust 
inputs. 


To facilitate the design and analysis of the automatic control systems. 

It is necessary to have the capability to describe numerically the interaction 
between structural, aerodynamic, and control forces. This interaction is 
described in terms of response and stability characteristics. This capability, 
in the form of a computer program system, has been assembled and packaged and 
is identified by the acronym I SAC, i.e. Interaction of Structures, Aerody- 
namics and Controls. The system is in reality an assembly of several 
programs tied together through a common data base. 

The following discussion is intended to explain the overall analysis 
capability of ISAC to a prospective user. The various mathematical, aerody- 
namic, structural, and control system theories underlying the program are 
given definitive treatment in the references. A brief statement of the 
related equations together with limitations on their applicability is given 
herein. 
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SYMBOLS 


A(s) 

B(s) 


coefficient of equations of motion which represent 
airframe motion 

undetermined coefficients of polynomial, equation 15 
coefficient of excitation 
coefficients of equation 15 


c 


[D] 



r^[f(oi)] 

Ha 

^66 

H(lo)) 


i 

K 

k 

M 

etc. 


reference chord length 

downwash, pressure differential influence matrix 

elements of D matrix 
exponential number 
Fourier transform of f(o)) 

servo transfer function 

surface servo transfer function 

servo-main structure transfer function 
frequency response function 

deflection of aerodynamic and structural elements 

/-T 

Kernel function 
reduced frequency, cw/ZU 
Mach number 

generalized mass matrices 
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g", o'" 


Si 


S 

T 


T^(s) 


U 

u 

w(x, y, 2 ) 

Wq 

X 


"i 


c 

X, y, z 
n, C, 6 
AP(x, y, z) 

6 

Y 

u 

COr 


number of zero crossings 

number of b coefficients 

reduced Laplace variable cs/2U 

matrix of generalized aerodynamic forces due to 
airplane motion and gust input 

2 

dynamic pressure, 1/2 pU 

ith modal displacement stress coefficients 

Laplace variable 

geometric transformation matrix 

transfer function relating airplane motion to 
control surface actuator valve input 

true airspeed 

excitation matrix 

downwash due to airplane motion 

downwash due to gust 
generalized coordinate matrix 

servo input command 

command to actuator valve 

Cartesian coordinates 

integration variable in equation 2 

pressure differential distribution 

control surface rotation generalized coordinate 

real part of Laplace variable 

circular frequency 

natural frequency of surface 
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to 


c 




a 

4>(x, y) 



natural frequency of structure 

airplane motion generalized coordinate 

RMS value, defined by equation 11 

shape functions from vibration analysis 
input power spectrum 

THEORETICAL BASIS 


In ISAC, the airplane is described in ter^s of a finite number of an 
infinite set of modes which represent the continuous structural/aerodynamic 
system. The ecuaticns of motion in terms of generalized coordinates, general- 
ized masses, springs, dampers and aerodynamic forces are essentially the same 
as those developed and discussed in reference 1 . The control systems a>^e 
treated as additional degrees of freedom. Control surface deflections are 
related to sensor inputs through transfer functions. 


Aerostructural Interface 


An analysis proceeds from modal characteristics obtained separately by 
any suitable vibration analysis program such as that contained in Nastran, 
reference 2, to the aero/structure interface program DLIN. In essense DLIN 
finds a transformation, [T] , utilizing the surface spline interpolation 
technique, reference 3, between deflections at structural nodes and deflections 
and slopes on aerodynamic boxes. 


h 



dh 

= [T] 

h 

dx 

a 



(I) 


Once available, the aerodynamic deflections and slopes are used in the 
computation of aerodynamic pressures. 


Generalized Aerodynamic Forces 


An unsteady aerodynamic lifting surface theory known as the Doublet 
Lattice method, reference (4), is used to compute the generalized aerodynamic 
forces and is contained in program DLAT. The downwash (velocity of air 
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particles normal to surface), at chosen reference points over the surface is 
equated to the surface integral of the differential pressure and an appropriate 
Kernel function is used to form the equation: 


w(x, y, z) 


[D] 


AP(x, y, z) 


where the elements of D are: 

*/ M) d^dn 

and the downwash in terms of sinsuoidal surface deflections is: 




dh 



w 

= U 

d 

+ io) 

h 

dx 

d 


( 2 ) 


The downwash arrays are computed internally in DLAT from the deflection 

A h 

and slope °”a arrays supplied by DLIN. Generalized aerodynamic forces 


dx 

matrices 


h 


a 


[Q(k, M)] = [4.]"^ [AP(k, M)] 


(3) 


are generated by DLAT for use in the equations of motion described in the 
following section. The program DLAT contained herein is based on the program 
described in reference 4, except that programing modifications such as those 
described in reference 5 have been made. These modifications take full advan- 
tage of the particular NOS operating system currently in use at LaRC to 
increase program efficiency from the standpoint of core/time usage. 


Equations of Motion 


The equations of motion employed to describe the vehicle are, as 
previously stated, similiar to those developed for small perturbation aero- 
el astic analyses in reference 1 . The reference axes shown in the sketch below 
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Reference Axes 


h(x, y, t) y)?(t) +^c|)^(x, y)6(t) 



are space fixed inertial and it is assumed that airplane motion relative to 
these axes is small so that the inertia tensor used to describe the system 
kinetic energy is unaltered. For clarity, the generalized coordinate matrix 
is partitioned into primary structure coordinates C , and control surface 

rotation coordinates 5 ^ and, in like manner the equations are separated into 

generalized forces in the primary structure modes lift and moment etc. and 
generalized control surface hinge moments. If it is assumed that the Laplace 
transform of the generalized aerodynamic forces, Q, exists, the equations of 
motion, without an automatic control system, are: 


M^^s^ + (1 + 


C + 


2 M 


q Q W- + H^^x^ 
G ^6 1 




s + 




6 = q Q^Wp + H;,~x. 
^“0 G 5o 1 


(4) 

(5) 


p 

The aerodynamic terms, Q on the right side of equation (4) and (5) represent 
gust inputs. The control surface hinge moments included on the right side of 
equation (5) result from an input command, x^ , to a hydraulic servo and a 

flexible surface, represented by the generalized force The generalized 

servo force acting on the structure is represented by the transfer function 
If a signal proportional to the primary structural motion and its 

derivatives is fed back to the servo as shown in the block diagram, figure 1, 
the servo input then consists of the sum of the feedback and a pilot command 
input 


X 


i 


= X. 



( 6 ) 


6 



where 




The matrices T^^p, Tp^ and <p represent sensor output to actuator input, 

sensor dynamics tranfer function and modal displacements and slopes respective- 
ly. When equation 6 is combined with equation 5 we get: 




? + 


Ml * 


= q QX + 

G 00 c 


(5a) 


Equation (5a) can be simplified by assuming an infinitely rigid surface and an 
irreversible servo. The generalized hinge moment becomes (surface rigid) 

and equation 5a becomes: 




\ 




6 = 


(5b) 


The substitution of equation (5b) into equation (4) yields: 


1^^s2 + (1 + iG)M^^c. 2 + qjM^ - 


55 


r -] 


^ 56 ^A 

^5 


q..QrWr. - + qj' 




?<S 


) 


«{««A 


(7) 


Equation 4 with 5a or 7 is used to describe the motion of the aircraft 
and the control system and can be written more concisely as follows: 


A(s) X = B(s)u 


( 8 ) 


where A(s) is the left hand side of either equations 4 and 5a or 7, i.e. a 
matrix of the coefficients of ^ and 6 which now are denoted and B(s) 

IS the matrix of the coefficients of the right hand side of 4 and 5a or 7, and 

77 contains W^. and x . 
u G c 

Equation 8 is essentially a set of simultaneous nonhomogeneous linear 
equations. The coefficients A(s) and B(s) consist of either rational 



polynomials of s of first or higher order or tabular functions of reduced 
frequency k and Mach number M. In either case, the primary tasks are to 
solve equation (8) for the generalized coordinates and to determine system 
stability. Whenever the closed form Laplace transform exists, standard pro- 
cedures are available for solving for the transfer functions 


X = A(s)"^ B(s) u 


( 9 ) 


and then obtaining the inverse X(t), and for determining the system stability 
from the roots of the characteristic equation |A(s) | =0. ISAC is designed to 
treat equation (8) both when the coefficients A(s) and B(s) are tabular 
functions of s=iu) or when they are functions of complex s. 

In the former case, solutions are obtained for J after first obtaining 
the frequency response functions 

H(i(jo) = Xiico) = A(io))'^ B(iw) (10) 

ujioj) 


Solutions found this way are valid only for stable systems (a result of the 
definition of the frequency response function), consequently it is essential 
to ascertain stability. Two methods of accomplishing this will be discussed 
later. 


If the forcing function, u, of equation (8) above is a nondeterministic , 
but Gaussian, random process, then two quantities, a (the RMS value) and N 

(the average number of positive going zero crossings), suffice to describe the 
response of the generalized coordinates to this excitation. These quantities 
are found from the frequency response function by the equations 


oo 


((jj)du) 


( 11 ) 


N 


0 





cj^|H(iu)) |^$^(u))du) 
0 


( 12 ) 


where $ (u) is the power spectrum of the excitation u . It is further 

U J 

assumed that gust inputs (atmospheric turbulence) can be represented by a one 
dimensional point spectrum and that all u- are uncorrelated. On the other 

J 

hand if the forcing function is deterministic, for example, a discrete gust or 
a step control input, a time history solution for x" is given by the inverse 
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Fourier transform, 


X(t) = F"^[H(ioj)u(iu)] 


(13) 


This computation is done in ISAC using Fast Fourier Transform techniques. 
Solutions for response quantities which are composed of a combination of the 
generalized coordinates and their derivatives are similarly obtained. As an 
illustration, consider the structural stresses that are generated during the 
airplane response to gust or control surface excitation. 

In ISAC, structural stresses are calculated by the modal displacement 
method 


a 


P 



(14) 


where is the stress at some station p or the stress in some structural 
P (i ) 

element and ^ is the component of the stress Op contributed by a unit 

displacement in the ith generalized coordinate X^. . X^ can be taken either 

from equation 10 or 13 depending on whether the excitation is random or dis- 
crete. To repeat, it is necessary that the system be stable for the frequency 
response function to be valid. Methods employed herein for determining system 
stability are discussed in the next section. 

As stated previously, active controls are employed for the alleviation of 
structural loads (stresses) which result both from atmospheric inputs (gust) 
and aircraft maneuvers. In ISAC, the equations of motion can describe only the 
small linear deviations from a rectilinear equilibrium flight path (gravity, 
thrust, drag forces are omitted) of an airplane compressed into essentially 
flat plate lifting surfaces. The stresses thus calculated are incremental 
loads which would have to be combined with the Ig loads for design purposes. 
Moreover, maneuver loads are calculated from time history response of the air- 
plane to a pilot command. The load factor imposed is indicated by the output 
of a sensor located at the airplane center of gravity. The incremental 
structural load (for example, wing root bending moment for a given maneuver 
load factor) is the product of the load factor, the ratio of the structural 
load, and the center of gravity load factor for the given control surface 
input. 


SYSTEM STABILITY 


System characteristic roots are identified by an eigenvalue iteration 
scheme STABCAR, unpublished, and by a determinant search routine PROOT. PROOT, 
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used herein for flutter analyses and root loci analyses is a digital program, 
presently unpublished, developed by the coauthor from the principles in 
reference 6. The many useful number sets obtained from the p-k solution are 

discussed in a later section. Both methods of finding the characteristic 
values yield equivalent results. Both are subject to a degree of unreliability 
in that convergence to desired roots is not assured. The p-k method, how- 
ever, will almost always provide answers if sufficient care is taken in 
providing good initial guesses. 

System stability can also be examined by determining the roots of the 
characteristic polynomial of equation (8), for which the aerodynamic forces 
are treated as functions of s. DLAT, the aerodynamic force generator within 
ISAC, provides aerodynamic forces only on the imaginary axis s = Kjj. 

However, it is feasible using the concept of analytic continuation to 
provide an approximation to the aerodynamic forces for points near the imagi- 
nary axis. Theoretical descriptions of this approach are presented in 
reference 7. In ISAC the s-plane (actually p-plane where p = cs ) approxi- 

'W 

mation for the aerodynamic forces is of the form: 

Ob 

Q(p) = + Ajp + A^p^ ^1+2'" 

(P + b. ) 


where the real matrices A^, . . • .AOj^ ^ 2 determined by solving for the 

elements of these matrices that best fit, in a least squares sense, curves 
through 


Q(p = ik^,Q(p = ik2), . . .,Q(p = ik^) 


1^'^b + 2 


The b are chosen to span the frequency range of interest. This form for 
the approximation is similar to that described in reference 7. 

A capability for transforming the aerodynamics from functions of k to 
functions of p is included in the Dynares Program. The generalized aerody- 
namic forces, Q(p), thus compiled are used in equation 8 in the same manner 
as the Q(k). Specifically, values of Q(p) for any arbitrary value of the 
argument s are found by evaluating equation 15. 


PROGRAM STRUCTURE 


ISAC relies heavily on the LaRC CDC computer and its particular operating 
system. Various programs are called and executed through NOS procedure files. 
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The basic structure of ISAC is shown in figure 2. The hub represents a data 
complex or specifically a random access file (TAPE9) from which and to which 
numerical data required by or generated by the various independent program 
modules is passed. The data complex can be accessed by the ISAC modules which 
read and/or store needed data. This capability allows convenient manipulation 
of and selection from a large amount of data. For example, the data base 
might contain several structural models and aerodynamic data at several Mach 
numbers for each of the structural models. The Dynares module can select from 
this data a particular structural model and a set of corresponding aerody- 
namics at a particular Mach number, thereby relieving the user of much of the 
effort of data storage and manipulation. Greater accuracy is achieved and 
fewer transcription errors are made in this way. 

Since a portion of the Data Complex Manager (DCM) is included in each of 
the modules, it is possible to execute each either separately, or with 
appropriate job cards, in a single sequential run. Generally, separate exe- 
cution is preferred so that inspection of intermediate results can be made. 
Accuracy judgments should be made on the intermediate results because aero- 
elastic analyses still involve a certain degree of engineering judgment in 
choosing which structural modes to include, which interpolation method yields 
the best fit, and how to panel the lifting surfaces for the aerodynamic 
calculations. 


NUMERICAL ANALYSIS 


Numerical analyses are carried out in three steps, (1) description of 
airplane configuration, (2) instructions for solution types required and 
(3) specification of results to be output. The configuration description 
includes structural paramete'^s, aerodynamic properties and control system 
definition. 

The structural parameters; generalized masses, generalized stiffnesses, 
damping ratios and modal frequencies are obtained from independent vibration 
analyses. These data may be stored on (TAPE9) using the DCM or read into 
Dynares on (TAPES). Dynares also has the option of computing generalized 
masses and stiffnesses from mode shapes and lumped mass data input through 
(TAPES). 

Generalized aerodynamic forces can be computed internally and stored on 
(TAPE9) or input directly to either (TAPE9) (DCM) or (TAPES) (DYNARES). 

DLIN, executed first, is used to calculate mode shape data for the aerodynamic 
boxes from the flexible structural mode shapes input on (TAPES). Wing 
geometry and paneling information is obtained from the same name list DLATINP 
(TAPE2) used by DLAT. The aerodynamic mode shapes are stored on (TAPE9) by 
selected DCM subroutines contained in DLIN. These data, identified as H14, H, 
and DH; are subsequently used in DLAT. Differential pressure distributions are 
calculated in DLAT for the flexible modal data either stored on (TAPE9) or read 
in on (TAPES). Rigid body modes, control modes and sinusoidal gust mode shapes 
are generated internally. Generalized aerodynamic forces defined by equation 
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(3) are obtained from the pressure distributions generated for a variety of 
reduced frequencies and Mach numbers specified in DLATINP (TAPE2). The 
generalized aerodynamic forces are stored on (TAPE9) for use in Dynares whereas 
the pressure distribution, i.e. Aerodynamic influence coefficients, (very 
large arrays) may be saved on magnetic tape for other uses if desired. 

After a description of the structural and aerodynamic properties are 
complied, the control system elements are input to Dynares through the name 
lists, INPUT, CONSYM, ACTINP, FILTIN, and SENLOC. Transfer functions 
representing sensors and filters which possess specific frequency response 
characteristics may be computed internally. In addition, sensors may be moved 
from one location to another by setting ISENSE=1 and supplying the sensor 
location in SENLOC namelist. 


To reiterate, the generalized aerodynamics forces are tabulated functions 
of Mach number and the reduced frequency k, and are based on a sinusoidal 
downwash. If, however, the downwash is nonsinusoidal , the aerodynamic forces 
are of questionable validity. To circumvent this difficulty, an approximation 
for Q(s = Y + i2Uk) deemed valid for nonnegligible y is achieved by trans- 


forming through a least squares fit the complex function Qp ^(k) + 




■“Real 

of the real argument k to a real function of the complex argument 
s = Y + 2UM and treating the result as if it were the Laplace transform of time 
c 

dependent generalized aerodynamic forces. This transformation is accomplished 
in the program by setting ISPLANE=1 and supplying values for the parameters 
NCOEF, BN, NPOLYC, ICOF. Following transformation of the aerodynamics, all 
solutions follow as before except that the p-k analysis is now a p-p 
analysis, i.e. as far as the user is concerned, stability calculations are 
performed in the same manner in either case. 


CONCLUDING REMARKS 


The program described herein represents a research tool for the 
aeroelastician who must include automatic controls systems in an airplane 
analysis or the controls engineer who needs to include the effects of flexible 
structures and unsteady aerodynamics in the design of automatic control systems. 
The program is limited in some respects for both parties. Structural stresses 
are not detailed enough for structural sizing nor are the details of the control 
system hardware faithfully represented. Nevertheless, knowledge of the im- 
portant interplay between the various parts of the total system can be provided 
through judicious parametric studies. 

As mentioned earlier, this program is tailored to the LaRC computing system 
and hence may be unusable directly on other computing systems. Nevertheless, 
the possibility exists that someone with comparable computer facilities may 
find use for ISAC. Consequently, this program is documented to relieve them 
from duplicating the effort described herein. 
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THE APPENDIX, ISAC INPUT DATA 


Included in the following tables are instructions for compiling input 
data decks and definitions of the parameters involved. Each program is treated 
in the chronological order generally used. However, as previously stated, the 
stand alone programs can be executed consecutively in a single run. Intra 
program communication media are defined in Figure Al. 

As presently constructed, the DLAT preprocessor, the DLIN, can accommodate 
noncoplanar horizontal lifting surfaces only. If vertical surfaces and bodies 
of revolution are to be included, the mode shape input to DLAT must be input 
through tapes and is not set up by DLIN. 

If NB^O, the body mode shapes must be input to DLAT directly. DLIN 
is not constructed to compute these data. The option for including additional 
modes will enable the user to input mode shape data for the body. This is 
accomplished in the following manner. NELAST, NR and NG are set equal to zero. 
Then NMD = NADD + NC and H14, H and DH are input for every box for all NADD 
modes. The NADD modes now will include all rigid body, elastic and gust modes. 
The control option, NC, will, however, still compute appropriate control surface 
mode shapes. The order of compiling card images input are given for each of 
the four ISAC programs in Table 1 through Table 4. Definitions of each of the 
input quantities are listed in tabular form following each deck listing. 

Results of a flutter analysis using ISAC are included in the appendix, 
table 5 and figures A2 through A4. The sample problem is based on the jet 
transport characteristics given in reference 1. The analysis can be reproduced 
by requesting the data files from the authors. 
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Data Complex Manager 
Table la 


PCM - Order Of Input . - 
File: INPUT 

The input file consists of a string of operation codes each followed by 
card images consisting of data required by that operation code as follows: 
Operation Codes: STORE, PRINT, WRITE, TOC, REWIND, END 

If Op Code is: 

a) TOC or REWIND, 

No further data required 
another OPERATION CODE. 

b) STORE, then add the following: 

TYPE 1 data 

TYPE 2 data (if needed) 

TYPE 4 data 

another OPERATION CODE 

c) PRINT or WRITE, then add the following: 

TYPE 1 data 

TYPE 2 data (if needed) 

TYPE 3 data (if needed) 
another OPERATION CODE 

The input file is terminated when the operation code, END is encountered. 

File: (TAPES) 

All formatted arrays to be stored on Data Complex in the order in which 
they are to be stored. This is determined by the user. 

File: (TAPE?) 

All binary arrays to be stored on Data Complex in the order in which they 
are to be stored. This also is determined by the user. 
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Definition Of Input . - 


This program performs data management and manipulation of Data Complex for 
program ISAC, built-in capabilities are as follows: 


1. 

STORE 

2. 

PRINT 

3. 

WRITE 

4. 

TOC 

5. 

REWIND 


Store data into a dataset in Data Complex from another 
dataset or from (TAPES), if formatted, or (TAPE7) if binary. 

Print out a specified array from Data Complex to file OUTPUT. 

Write out a specified array onto (TAPES) if formatted, or 
(TAPES) if binary. 

Print out Table of Contents of entire Data Complex. Or TOC 
N print out Table of Contents of dataset N. 

N Rewind (TAPE N). 


The array input and output can be either formatted or binary. If formatted, 
the format must also be specified on file INPUT. 

STORE, PRINT, and WRITE, Operations require additional input - as follows 

(TOC and REWIND do not): 

Type 1: CODENAME, INPUT PARAM, OUTPUT PARAM, NWORDS, NRECORDS 


END 


Type 2: Format for reading from (TAPES) if needed. 

Type 3: Format for output-PRINT or WRITE - if needed. 

Type 4. Header for storing - if needed. 

Definitions: 

CODENAME Code name of array to be processed as follows: 


GMASS 

(Generalized Masses) 

FREQ 

(Natural Frequencies) 

AERO 

(Red. Freq + Aero Table) 

SPLANE 

(S-Plane Coefficients) 

LOADS 

(Loads-Shear , Moment, Torque) 

SENDEE 

(Sensor Deflections) 

DAMPINGS 

(Structural Dampings) 

DISP 

(Modal Displacements) 

MODES 

(Selected Mode Shapes) 

HHD 

(H14, H34, DHl Arrays) 

SPLINE 

(SPLINE COEFFICIENTS) 

AIC 

(AIC MATRIX) 
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INPUT PARAM A number indicating from where and how array is to be obtained. 


OUTPUT PARAM 

NWORDS 

NRECORDS 


=N>0 Read from data-set N. 

=-l Read from (TAPES) in formatted form. 

=-2 Read from (TAPES) in binary form. 

A number indicating where and how array is to be output. 

=N>0 Store on data-set N. 

=-l Output (PRINT, WRITE, or PUNCH) in formatted form. 

=-2 Output (PRINT, WRITE, or PUNCH) in binary form. 

Number of words in one record of array to be processed 
(input only if array is not already on Data Complex) 

Number of records used to comprise entire array (input only 
if array is not already on Data Complex) 
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DCM-Sample Problem . - 


Suppose it is desired to determine the contents of a dataset N of 
a data file XXX, then modify a portion of the set N update XXX and finally 
check the content. The following control deck and input file would be used. 

JOB ... 

USER ... 

CHARGE ... 

MAP, OFF. 

GET,DCM/UN=887010C. 

DCM. 

EDIT, (TAPES). 

RENAME, (TAPES )= (TAPES). 

DCM 

REPLACE, (TAPE9)=XXX. 

7/8/9 

TOC 

WRITE 

GMASS N -1 

END 

(4E15.7) 

END 

7/8/9 

...EDIT COMMANDS... 

7/8/9 

STORE 

GMASS -1 N 

END 

(4E15.7) 

MODIFIED MASSES 

TOC N 

END 

6/7/8Z9 
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Aerostructural Interface 


Table 2a 


DLIN - Order Of Input . - 
File: TAPE3 

Starting 



Description 

Format 

Column 

Data 1 

HEADER 

8A10 

Column 1 

Data 2 

NVIB 

Free field 

Column 1 

Data 3 

DLINPT Namelist 

Name list 


Data 4 

MODEPLT, Namelists (REPEAT FOR EACH SET OF MODE SHAPES) 

Data S 

DCM input-output parameters 

Free field 

Column 1 

Data 6 

1=1, NSECTNS 



(NSECTNS CARDS) 

ISS(1,I),ISS(2,I),IPLATE(I), 

Free field 

Column 1 


X0(I),Y0(I),R0(I) 



Data 7 

1=1, NSECTNA 



(NSECTNA CARDS) 

IAS(1,I),IAS(2,I),NSS(I),ID(I) 

Free field 

Column 1 

Data 8 

1=1, NNODES 



(NNODES CARDS) 

TAB(I,l),TAB(I,2),Nji)DE(l), 

Free field 

Column 1 


IDF(1,I),IDF(2,I) 

File: TAPEl 

Data 1 IM=1,NELAST 

SKIP NSKIPL {Lines if formatted, records if binary) 
I=l,NNODET 

READ(DEFL(1,J),J=1,NDF) if formatted (18X,3E18.6/18X,3E18.6) 
This is unnecessary if these deflections have previously been put into the 
Data-Complex (TAPE9) under the Codename - DISP. 

File: TAPES 

All TAPES data for DLAT must 
also be input to DLIN. (See 
input for DLAT), 
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Table 2b 


Definitions of DLIN input data 

FILE: TAPE3 .- 

Data 1 

HEADER - any user message 


Data 2 

NVIB - the NVIB parameter tells the program how the deflection data 
are coded on TAPEl. 

If NVIB=NASTRAN, the program defaults NSKIPL=7 and reads in the 
deflection data according to the format (18x,3E18.6/18x,3E18.6) 

If NVIB=SPAR, the program defaults NSKIPL=0 and reads in the 
deflection data as an unformatted binary form. 

If NSKIPL IS specified in DLINPT, that value will overide the above 
program defaults. Be careful to maintain consistency. 

If NVIB IS anything other than NASTRAN or SPAR, NSKIPL will be the 
value specified in DLINPT and the read format must be given. For example, 
two cards are required for NVIB: 

ANYTHING Card 1 

(2E15.6) Card 2 
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DATA 3, DEFINITIONS FOR NAMELIST /DLINPT/ 

NAMELIST /DLINPT/ NSKIPL,NSECTINS,NSECTNA,NDF,NNODES,NNODET 
C , ISYM , MOD E , NPMX , I PLTBOX , SCX , I PLTNOD , I PR I NT , I PL IMS 
C .INTER, IDCM,ITOC,IBAUD 


NSKIPL 


NSECTNS 

NSECTNA 

NDF 

NNODES 

NNODET 

ISYM 

M0DE(I)=1 


NPMX 

IPLTBOX=l 

=0 

sex 


SCY 

NOTE: 


IPLTN0D=1 

IPRINT=0 
= 1 

IPLTMS=0 
= 1 
=-l 


NO. OF LINES (OR RECORDS) TO BE SKIPPED AT BEGINNING 
OF EACH MODAL SET OF INPUT DATA ON TAPEl. 

(DEFAULT FOR SPAR = 0, FOR NASTRAN = 7) 

NUMBER OF STRUCTURAL SECTIONS 

NUMBER OF AERODYNAMIC SECTIONS (OR PANELS) 

NUMBER OF DEGREES OF FREEDOM OBTAINED FROM VIBRATION 
PROGRAM 

NUMBER OF NODES TO BE SELECTED FROM TOTAL OBTAINED 
FROM VIBRATION PROGRAM 

TOTAL NUMBER OF NODES OBTAINED IN VIBRATION PROGRAM 
MODE SHAPE SYMMETRY CODE FOR SPLINE COEFFICIENTS 
(IF 0, UNSYMMETRIC, IF 1, SYMMETRIC, IF -1, ASYMMETRIC) 
INDICATES THAT THAT MODE WILL BE SELECTED FROM 
IN DOING NODAL SELECTIONS IN 0VERLAY(2 ,0)-N0DESEL 
(DEFAULT = 1) 

MAXIMUM NUMBER OF NODES IN ANY ONE STRUCTURAL 
SECTION 

PLOT AERODYNAMIC BOX LAYOUT 

DO NOT PLOT BOX LAYOUT 

SCALE FACTOR FOR X-VALUES IN BOX LAYOUT 

=XMAX/FSX, WHERE XMAX IS AN APPROXIMATE 

MAXIMUM X VALUE AND FSX IS THE FRAME SIZE IN 

THE X DIRECTION. 

SCALE FACTOR FOR Y-VALUES IN BOX LAYOUT 
=YMAX/FSY 

IN ORDER TO HAVE NO DISTORTION DUE TO DIFFERENT 
SCALE FACTORS FOR X AND Y, SET SCX=SCY= 
MAX(XMAX/FSX,YMAX/FSY) 

(DEFAULT VALUES FOR SCX AND SCY=INT(MAX(X,Y)/10. ) 

PLOT STRUCTURAL NODE POINTS WILL BE PUT ON 

BOX LAYOUT IF IPLTB0X=1,0R PLOTTED ALONE IF IPLTBOX 

=0. sex AND SCY WILL STILL USED. 

PRINT ORIGINAL MODE SHAPES ON FILE OUTPUT 
(DEFAULT) 

PRINT ORIGINAL MODE SHAPES ON FILE TAPE6. 

MODAL DATA WILL BE COMPUTED FOR AERO. BOXES 
(DEFAULT) 

MODAL DATA FOR AERO. BOXES IS COMPUTED; AND 
MODAL DATA FOR PLOT POINTS IS COMPUTED AND PLOTTED. 
COMPUTES AND PLOTS MODAL DATA FOR PLOT POINTS 
USING SELECTED DATA FROM TAPE9 (DATA COMPLEX) IF 
AVAILABLE FROM A PREVIOUS RUN. 



INTER 

=0 

BATCH TYPE RUN. ADDITIONAL MODEPLT DATA 
READ IN FROM TAPES. 

(DEFAULT) 


=1 

INTERACTIVE TYPE RUN, ADDITIONAL MODEPLT DATA 
READ IN FROM FILE INPUT. 

IDCM 

=1 

ALLOWS USER TO ACCESS DATA-COMPLEX 

ITOC 

=1 

PRINT TABLE OF CONTENTS FOR DATA-COMPLEX 

I BAUD 


BAUD RATE (DIVIDED BY 10) TO TERMINAL WHEN TEKTRONIX 
PLOT ROUTINES ARE BEING USED 

END 

NAMELIST DLINPT 
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o o o 


DATA 4, DEFINITIONS OF VARIABLES IN NAMELIST/ MODEPLT/ 

(IF IPLTMS = +1 OR -1, MODEPLT NAMELISTS MUST BE SUPPLIED. 

REPEAT FOR EACH SET OF MODESHAPE PLOTS.) 

NAMELIST /MODEPLT/ NSURF,XLE,YLE,NLEP.XTE,NTRP.NSPN 
, NCHD , NSSP , SC XMS , SC YMS , SCZMS . SC ZD EF , XV , YV , ZV 
,ALPHA,BETA,XTRANS,NE,MODEMS,XOR,IPLAN,IBOTH 
,INODELN,ISTOP 

NSURF NUMBER OF SURFACES TO BE PLOTTED 

XLE(I,J) X,Y - COORDS OF I"TH LEADING EDGE POINT ON J"TH 
SURFACE 

YLE(I,J) NEEDED TO DEFINE LEADING EDGE (OF SURFACE NUMBER J) 

NLEP(J) NUMBER OF LEADING EDGE POINTS NEEDED TO DEFINE J"TH 

SURFACE 

XTE.YTE DEFINED THE SAME AS ABOVE EXCEPT FOR TRAILING EDGE 
NTRP(5) NUMBER OF TRAILING EDGE POINTS NEEDED TO DEFINE J"TH 
SURFACE 

NSPN(J) NUMBER OF SPANWISE POINTS TO BE PLOTTED FOR J"TH 
SURFACE 

NCHD(J) NUMBER OF CHORDWISE PNTS TO BE PLOTTED FOR J"TH 
SURFACE 

NSSP(J) STRUCTURAL SECTION NUMBER WHICH CORRESPONDS TO J"TH 
SURFACE 

SCXMS SCALE FACTOR IN X-DIRECTION OF VIEWING IMAGE 

SCYMS SCALE FACTOR IN Y-DIRECTION OF VIEWING IMAGE 

SCZMS SCALE FACTOR IN Z-DIRECTION OF VIEWING IMAGE 

SCZDEF SCALE FACTOR FOR DEFLECTIONS, USED TO INCREASE OR 

DECREASE AMOUNT OF VIEWED DEFLECTION 
XV,YV,ZV COORDS OF VIEW POINT FOR PROJECTIVE PLOTS WITH 
RESPECT TO IMAGE SIZE. 

ALPHA ANGLE OF ROTATION ABOUT Z-AXIS OF VIEWING IMAGE 

BETA ANGLE OF ROTATION ABOUT Y-AXIS OF VIEWING IMAGE 

XTRANS AMOUNT OF X-TRANSLATION OF VIEWING IMAGE 

NE NUMBER OF TOTAL ELASTIC MODES TO BE PLOTTED 


MODEMS(I)=l PLOT MODE SHAPE FOR ELASTIC MODE I 

=0 DO NOT PLOT MODE SHAPE FOR ELASTIC MODE I 
XOR,YOR LOCATION OF PLOT FRAME ORIGIN, USED TO SHIFT 
PLOTTED IMAGE IN PLOT FRAME 

IPLAN =1 PLOT UNDEFLECTED PLANFORM AND CONNECTING VECTORS 
TO DEFLECTED PLANFORM. 

=-l PLOT UNDEFLECTED PLANFORM WITHOUT CONNECTING VECTORS 



=0 

DO NOT PLOT PLANFORM. 



(DEFAULT=1) 

I BOTH 

= 1 

PLOT BOTH SIDES OF VIEWING IMAGE 


=0 

DO NOT PLOT BOTH SIDES 

INODELN 

= 1 

PLOT NODE LINE ALONG WITH DEFECTED PLANFORM. 


=-l 

PLOT NODE LINE ONLY. 


=0 

DO NOT PLOT NODE LINE. 


(DEFAULT = 1) 
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ISTOP =0 CONTINUE READING IN /MODEPLT/ NAMELISTS AND 
PLOTTING MODESHAPES 

=1 STOP READING IN /MODEPLT/ NAMELISTS. THIS MUST 
BE SET = 1 IN LAST /MODEPLT/ TO BE READ IN. 

END NAMELIST MODEPLT 


Data 5 

DCM Input-Output Parameters are card images of the following type: 

CODENAME, INP,OTP 
where 

CODENAME: DISP, MODES, HHD , SPLINE, or END 

DISP Modal displacements (deflections) from vibration program. 

MODES Selected MODE shapes (set up) by overlay Nj^DESEL. 

HHD Quarter chord deflections and slopes set up by overlay DEFLEC. 

SPLINE SPLINE COEFFICIENTS, used for interpolation, set up by overlay 

INTERP. 

END Ends DCM input. 

INP An INTEGER which indicates from where its corresponding 

array is to be obtained. 

=-l read from TAPEl input file 
= 0 not used as input to program 

=n>0 read in from data-set n of Data Complex. 

OTP An INTEGER which indicates on which data-set of the Data 

Complex the array is to be stored. 

=0 Not stored 

=n>0 Stored on data-set n. 

Note: For DISP, INP can be -1 or n>0 

OTP can be 0 or >0 

For MODES, INP = OTP = n>0 

For HHD, INP = OTP = n>0 

For SPLINE, if there are any PLATE type modes, then INP=0TP=n>0; 
otherwise INP=0TP=0. 



DATA 6, STRUCTURAL MODEL DESCRIPTION 


ISS(1,I) BEGINNING NODE NUMBER OF STRUCTURAL SECTION I 

ISS(2,I) ENDING NODE NUMBER OF STRUCTURAL SECTION I 

IPLATE(I)=1 INDICATES THAT STRUCTURAL SECTION I IS PLATE TYPE 

IPLATE(I)=0 INDICATES THAT STRUCTURAL SECTION I IS BEAM TYPE 

XO(I) X COORD OF TRANSLATED AXIS OF SECTION I 

IF IPLATE(I)=1 - (XTE+XLE)/2 FOR THAT SECTION 
YO(I) Y COORD OF TRANSLATED AXIS OF SECTION I 

IF IPLATE(I)=1 - YROOT FOR THAT SECTION 
RO(I) ROTATION ANGLE OF ELASTIC AXIS OF SECTION I 

IF IPLATE(I)=1 - 1/B WHERE B=CROOT(I)/2 

DATA 7, AERODYNAMIC MODEL DESCRIPTION 

IAS(1,I) BEGINNING BOX NUMBER OF I-TH SECTION 
IAS(2,I) ENDING BOX NUMBER OF I-TH SECTION 

NSS(I) STRUCTURAL SECTION NUMBER WHICH CORRESPONDS TO I-TH 

AERODYNAMIC SECTION 

ID(I) =XY, IF I-TH SECTION IS HORIZONTAL 

=XZ, IF I-TH SECTION IS VERTICAL 

END TAPES DEFINITIONS 

DATA 8-. Structural nodes (see Sketch Al) 

TAB = COORDINATES OF NODES 

NODE(IN) = NODE number to be selected 

IDF(1,IN) » First Degree of Freedom to select (SEE NOTE) 

IDF(2,IN)= (Optional) Second Degree of Freedom to select (SEE NOTE) 



SKETCH Al NOMENCLATURE STRUCTURAL SECTIONS 
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FILE: TAPEl." Flexible mode shapes 

DEFL(I,J) Is the J'th deflection for the I'th structural NODE 

NOTE: THE SHAPE FUNCTIONS DEFL ( I, J), OUTPUT FROM NASTRAN CONTAIN SIX DEGREES 

OF FREEDOM FOR EACH NODE. THE IDF PARAMETER SPECIFIES WHICH TWO OF 
SIX ARE TO BE USED. 

DLIN - Sample Problem. - The following job illustrates how DLIN might be 
executed separately to store aerodynamic downwash data on TAPE9 and/or obtain 
a structural NODE plot. 

DLINB, T200 
USER,XXXXXXX,XXXXXXX. 

CHARGE, XXXXXX,LRC. 

MAP, OFF. 

GET,TAPE9=BIS9. 

GET,TAPE1=NASBISC. 

GET,TAPE3=DLIBC. 

GET,TAPE5=DLABC. 

ATTACH , FTNML I B , LRCGOSF/UN=L IBRARY . 

GET,ACP0LIB/UN=887010C. 

GET,LG0=DLINBN/UN=887010C. 

LD SET ( L I B=ACPOL I B/ FTNML I B/ LRCGOSF ) 

LGO. 

PL0T.VARIAN(FSH=11. ,FSV=8.5,AUT0(0)) 

REPLACE, TAPE9=BIS9. 

DAYFILE,DLINDAY. 

REPLACE, DLINDAY. 

EXIT. 

DAYFILE, DLINDAY. 

REPLACE, DLINDAY. 

END OF FILE 

(The Data Files NASBISC, DLIBC and DLABC have been previously stored, DLIBC 
and DLABC could be card input, using COPYL, INPUT, XXX) 

ACPOLIB IS a library of data management routines and plot routines that may be 
required in various modules of ISAC. 


BINXX XXXXX 

Omit if B1S9 has not been created 


Generalized Aerodynamic Forces 
Table 3a 

PLAT - Order Of Input . - 

File: (User Supplied Name) will be accessed as TAPES by DLAT 
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Description 


Format 


Starting 

Column 


Data 1 

HEADER 

8A10 


Column 2 

Data 2 

Namelist $DLATINP$ 

Namel ist 

Column 2 

Data 3 

DCM input-output parameters 

Free field 

Column 1 

(Reads MODE shapes from 

TAPE9 if IDCM=1, or from TAPED, 

Ser 1 if IDCM=0) 

Data 4 

1=1, NP 





(XCAP(J),J=1,4),(YCAP(J),J=1,2) 

Free field 

Column 1 


(ZCAP ( J ) , J=1 ,2 ) ,NSPAN ,NCH0RD , 





COEF 

Free field 

Column 1 


(TH(I),I=1,NCH0RD) 

Free field 

Column 1 


(TAU(I),I=1,NSPAN) 

Free field 

Column 1 

Data 5 (if NB/o) 

1 = 1, NB 





ZSC,YSC,NF,NZ,NY,COEFF, 





MRK(1),MRK{2) 

Free field 

Column 1 


(F(I),I=1,NF) 

Free field 

Column 1 


(RAD (I), 1=1, NF) 

Free field 

Column 1 

Data 6 

(LIM(I,1),LIM(I,2),I=1,NSTRIP) 

Free field 

Column 1 

(NADD additional MODE data can be read in here instead 

of or in 

1 addition to 

that on TAPED IF NADD = 

= NMOORS - (NELAST + NR + NC)=0. 

NADD IS COMPUTED 

INTERNALLY. 





Data 7 (if NADD/O) 






1 = 1, NADD 





(H14(J),J=l,NBOX) 

(6F10.0) 

Column 1 


(H(J), 0=1, NBOX) 

(6F10.0) 

Column 1 


(DHl(J), 0=1, NBOX) 

(6F10.0) 

Column 1 

Data 8 (if NC/O 

1=1, NC 





NBC,NCST,SWHL,ICTYPE 

Free field 

Column 1 


(ICB0X(IK),IK=1,NBC) 

Free field 

Column 1 

Data 9 (if NBE/O) 

1=1, NMODES (include all MODES shapes)"] 




BQ(I),I=l,NBOXB 


REPEAT FOR 


H(I),I=l,NBOXB 


EACH 

OF 


DH1(I),I=1,NB0XB 


NBE BODIES 


DH2(I),I=l,NBOXB 





Table 3b 


Definitions For PLAT Input . - 

Namelist - /DLATINP/FMACH,ACAP,REFCHD,NDELT,NP,NB,NRF,RFREQ,NSTRIP,NSV,NBV, 

NMODES , REFSPN , NBOX , NSNC , XZERO , XCG , I PR I NT , NR , NELAST , NC , IDCM , ITOC , 
IRIG 

FMACH Mach number 

ACAP Reference area 

REFCHD Reference length K=(REFCHD*0MEGA)/(2.*U) 

NDELT Symmetry Flag (=1 Symmetric ,=-1 Antisym,=0 No-Sym) 


27 



NP 

NB 

NRF 

RFREQ 

NSTRIP 

NSV 

NBV 

NMODES 

NDATA 

REFSPN 

NBOX 

NSNC 

XZERO 

XCG 

IPRINT 

NR 

NELAST 

NC 

IDCM 
I TOC 
IRIG 


DATA 3 

DCM input-Output Parameters are card images of the following type: 

CODENAME, INP,I0UT, (AICTAPE) 
where 

CODENAME Is either HHD,AER0,AIC, or END 

INP=n>0 Is an integer which indicates from which data-set in Data- 
Complex array is to be read. 

=0 Not read from Data-Complex 

I0UT=n>0 Is an integer which indicates into which data-set array is 
to be stored. 

=0 Not stored on Data-Complex. 

AICTAPE Filename of optional file to store AIC matrix on, other than 
TAPE 9. 

NOTE: for HHD, INP=n^0,I0UT=0 

for AERO, I0UT=n^0, INP= 

for AIC, if AICTAPE is blank, and either INP or lOUTj^O, then TAPE 9 
will be used for storage. 

DATA 4 

(XCAP(J),J=1,4) The panel edge coordinates 

(YCAP(J),J=1,2) 

(JCAP(J),J=1,2) See Sketch 2. 


Total number of panels 
Total number of bodies, see text 
Number of reduced frequencies 
Array of reduced frequencies 
Total number of strips 

Total number of strips on vertical panels on Y=0 
Total number of boxes on vertical panels 
Total number of MODES (elastic+rigid+control+NADD) 
Print flag for slender body solution 
= 1.0 

Total number of boxes on lifting surface panel 

Max NSPAN*NCH0RD value for any panel 

Gust reference station 

Rotation point for rigid body pitch MODE 

Print flag for B0,W and CP matrices 

Number of rigid body MODES 

Number of elastic MODES to be read from DATA-COMPLEX 
Number of control MODES shapes computed internally 
Allows user to access DATA-COMPLEX 
Print table of Contents for DATA-COMPLEX 
Array for rigid body Modes (IRIG (I), 1=1, NR) 

(1) = Z-Displacement (PLUNGE) 

(2) = 0-Rotation (PITCH) 

(3) = Z-Velocity W 

(4) = Angular Velocity Q 
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XCAP(1),YCAP(1),ZCAP(1) 

XCAP(2),YCAP(1),ZCAP(1) 

XCAP(3),YCAP(2),ZCAP(2) 


Inboard leading edge 
Inboard trailing edge 
Outboard leading edge 




(X^.Yj.Z 


0 

T 

(x,y,z) 

NSPAN 

NCHORD 

COEF 

THETA 

TAU 



SKETCH A2. AERODYNAMIC PANEL GEOMETRY 

corresponds to TH input data, Chord Ratio 
corresponds to TAU input data. Span Ratio 
corresponds to (XCAP,YCAP,ZCAP) input data. 

Number of spanwise divisions for panel considered. 

Number of chordwise divisions for panel considered. 

Scale factor for panel deflection modes. If no scale 
factor is desired, must be set to 1.0. 

Chordwise divisions in fraction of chord. Usually varies 
from 0.0 at leading edge to 1.0 at trailing edge. 

Spanwise divisions in fraction of panel span. Usually 0.0 
at inboard chord and 1.0 at outboard chord. 
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DATA 5 


Z-coordinate of body-axis. 
Y-coordinate of body axis. 
Number of divisions on body. 


ZSC-Z 

c 

YSC-Y 

c 

NF 

NZ If "1", the program will allow for doublets in the Z- 

di recti on and thus will be able to account for Z-or 
upwash. If not, set NZ = 0. 

NY If "1", program will allow for Y-doublets and thus Y-or 

sidewash. 

COEF Scale factor for body deflection mode. 

MRK Each body has associated with it interference lifting 

surfaces. MRK(l) represents the first box on these 
surfaces, and MRK(2) is the last box on the body 
interference lifting surface panels. Interference 
lifting surface panels are input after the regular 
lifting surfaces are input. 

F Endpoints of body divisions (starts at leading edge and 

ends at trailing edge). F=x coordinates of body 
element end points. 



Sketch A3. Body Geometry 
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DATA 6 
LIM 


DATA 7 
NADD 


H14,H 

DH1,DHZ 


DATA 8 

NBC 

NBE 

NBOXB 

NCST 

SWHL 

ICTYPE = 0 
= 1 
= 2 

ICBOX 
DATA 9 
BQ(I) 

H(I) 

DH1(I),DH2(I) 


Summation limits for chordwise integration. Usually LIM 1 
is the first box in a strip and LIM 2 is the last box where 
all boxes are numbered consecutively. 


Number of additional modes, NADD = NMODES-(NELAST + NR 
+ NC) NADD is computed by program, not input. The 
additional modes can be any arbitrary mode shapes. If all 
NELAST mode shape are to be input via TAPE 5 set NELAST 
= 0 . 

Deflections at 1/4 chord and 3/4 chord of box middle. 

1st and 2nd derivative with respect to streamwise direction. 
See Reference. 


Number of boxes on control surface. 

Number of slender bodies. 

Number of boxes on each slender body. 

Number of strips on each control surface. 

Sweep angle of hinge line or X-location of hinge, depending 
on the type of control surface. 

For trailing edge control. (Default) (Input Sweep Angle). 
For leading edge control, (input Sweep Angle). 

For all movable surface control. (Input hinge location). 

Box numbers of boxes on control surface, in ascending order. 


Integration matrix elements, i.e. area of ith box times 
deflection at quarter chord. 

Deflection at three-quarter chord of the ith box. 

First and second derivative of deflection shape at the 
three-quarter chord of the ith box. 


PLAT - Sample Input .- The Doublet Lattice can also be executely independently 
provided downwasti data has been previously stored or is input on TAPE 5. 
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Control Deck 


DLI^TD,TSQ9. BIHKX XKXXX 

UGER,XXXXXXX,X}IXXXXX. 

CP^RUE,XXXXXX,LRC. 

MAP, OFF. 

G£T,TAPE9«BIS5. 

GFT.TAPSSaDLftBC. 

G£T,LG0>DLATB»l/Un«8878ieC. 
ftTTftCHf FTHHLIB/UH-LI8RARY. 

GET , fiCPOL I B ^UN=8870 1 0C. 

GET,DYLIC/UN=387010C. 

LDSET(LI8='DYLIB/ACP0LIB/FTN«LIB. 

LGC. 

REPLACE, TAPE9*BIS9. 

DAYFILE.DLATDAY. 

REPLACE, DLATDAY. 

EXIT. 

DAYFILE,DLATOAY. 

REPLACE, DLATDAY. 

EKD OF FILE 


STABILITY AND RESPONSE ANALYSES 


The following table describes the procedure for executing the program 
DYNARES as a stand alone program. DYNARES will be utilized to a far greater 
extent that will either DLIN or DLAT since the latter programs generate 
configuration characteristics which may only be called upon once for the 
geometry and mode shape data and once each Mach number desired (DLAT only) 
whereas DYNARES will be used to perform stability and response analyses for a 
large number of variations of flight conditions and control system parameters. 


Table 4a 


Dynares, order of TAPE2 input data 


Description Format 


Data 

1 


NOCASES 

Free field 

Da ta 

2 


HEADER 

8A10 

Data 

3 


$ INPUT. . .$ 

Namelist 

Data 

4 

(If IDCM=1) 

DCM parameters 

Free field 

Data 

5 

(If ISENSE=1) 

$ SENLOC. . .$ 

Namelist 

Data 

6 

(If ICSPREAD/0) 

(CS(J,I),J=l,NM-Nc),I=NS) 

Free field 

Data 

7 

(If IC0NSYS=1) 
and if KASE*>1, 
ICHANG1=1) 

$ CONSYM. . .$ 

Namelist 
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Data 8 (If IC0NSYS=1 
and if KASE>1, 
ICHANG2=1) 

1=1, NC 
J=l, NS 
$ FILTIN. . 

.$ 

Namelist 

Data 9 (If IC0NSYS=1 
and if KASE>1, 
ICHAN6E3=1) 

$ ACTINP. . 

.$ 

Namelist 

Data 10 (If ISELECT=1) 

$ SELECT. . 

.$ 

Namel i st 

Data 11 (If any PLOT 

options are on) 

$ PLTSEL. . 

.$ 

Namelist 


If in Batch mode, repeat Data 2 - Data 10 for each case. In the 
Interactive mode, (NOCASES = -1), Data 3 through Data 11 are input via the 
terminal after the first case. Furthermore, the user will be asked at the 
beginning of each case whether he/she wishes to STOP. Answer 1, for yes. 
Depress Carriage Return) 0 (or CR) for no. 

In both modes, after the first case, only changes to data need be input. 

*KASE is a program computed index for tracking the case number. 
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Table 4b 


Dynares - definition of TAPE2 data. 

EXPLANATION OF NOCASES 

NOCASES NUMBER OF CASES TO BE RUN DURING EXECUTION OF JOB 

IF < 0, THE INTERACTIVE VERSION OF THE PROGRAM 
WILL BE IN EFFECT AND AFTER FIRST CASE IS RUN 
ALL ADDITIONAL INPUT WILL BE FROM TERMINAL 

EXPLANATION OF PARAMETERS IN NAMELIST INPUT 

NOTE - I PARAMETERS (LIKE IGUST,IFLUT,ETC) ARE ALL ON-OFF OPTIONS 

1 - INDICATES ON 
0 - INDICATES OFF (DEFAULT) 

N PARAMETERS INDICATE "NUMBER OF" 

IDCM=1 ALLOW USER TO MANAGE AND MANIPULATE 

=0 RETAIN AUTOMATIC PROGRAM MANAGEMENT 

IT0C=1 PRINT OUT TABLE OF CONTENTS OF DATA 


DATA 3, DEFINITIONS OF VARIABLES IN NAMELIST/ INPUT/ 

NAMEL I ST/ INPUT/ NM , NR , NC , NK , C , XCG , NPMX , NSECTNS , NNODES , 
ISS,IPLATE,XO,YO,RO,IMAT,NMASS,KFIT,ISPLANE,NCOEF,BN, 

NPOLYC , ICOF, I SENSE , NS , INTERP , I SYM , NGV , ELVN , I VANE , 

ISCREAD , ISELECT, ICHANGl , ICHANG2 , ICHANG3 , ICSACT , IFORCE , 
NFORCE,IFLUT,KVAR,NV,DV,VO,HO,RHOO,NIT,NCUT,NFINE, 
EPSI,I0PT1,I0PT2,ITRACE,P1I,P2I,IPS,DEN,I0PT,XMACH, 
C0NFAC1,C0NFAC2,C0NFAC3,IGAIN,DELGAIN,ISTABCR,IZER0ES, 
SS,OMR,IGUST.NG,NGR,RHO,NINT,RFLOW,RFCUT,GL,UU, 
NVEL,IVSLF,ISOUT,NLOAD.ITHRF,MINT,TMAX,IPR,KFVN, 
NGTH,IFRQOUT,NCP,IPRINT,TFAC,IAPLT,ISPLT,IPKPLT,ITHPLT. 

EXPLANATION OF PARAMETERS IN NAMELIST INPUT 


DATA COMPLEX 
OF DATA COMPLEX 
COMPLEX 


Definition of Structural and Aerodynamic model 

NM 
NR 


NUMBER OF MODES(RIGID .ELASTIC, AND CONTROL) 
NUMBER OF RIGID BODY MODES 
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NC 

NK 

C 

XCG 

NPMX 

NSECTNS 

NNODES 
ISS(1,I) 
ISS(2,I) 
IPLATE(I) =1 
IPLATE(I) =0 
XO(I) 

YO(I) 

RO(I) 

IMAT =1 

=0 

NMASS 

KFIT 


ISPLANE=1 


=2 


NCOEF 

BN(I) 

NPOLYC(J) 

IC0F(N,J)=1 

=0 


ISENSE 

NS 


NUMBER OF CONTROL MODES 
NUMBER OF REDUCED FREQUENCIES 
REFERENCE CHORD 

X-COORD OF THE ROTATION POINT OF THE RIGID BODY 
PITCH MODE 

MAXIMUM NUMBER OF NODES IN ANY ONE STRUCTURAL 
SECTION 

NUMBER OF STRUCTURAL SECTIONS 

DATA THE NUMBER OF NODES AND BOXES ARE DIVIDED INTO 
TOTAL NUMBER OF NODES IN ALL STRUCTURAL SECTIONS 
BEGINNING NODE NUMBER OF STRUCTURAL SECTION I 
ENDING NODE NUMBER OF STRUCTURAL SECTION I 
INDICATES THAT STRUCTURAL SECTION I IS PLATE TYPE 
INDICATES THAT STRUCTURAL SECTION I IS BEAM TYPE 
X COORD OF TRANSLATED AXIS OF SECTION I 
IF PLATE(I)»1,» (XTE+LE)/2 FOR THAT SECTION 
Y COORD OF TRANSLATED AXIS OF SECTION I 
IF IPLATE(I)*1,» YROOT FOR THAT SECTION 
ROTATION ANGLE OF ELASTIC AXIS OF SECTION I 
IF IPLATE(I)»1,= 1/B WHERE B=CR00T(I)/2 
INPUT LUMPED MASSES AND INERTIAS FROM WHICH 
GENERALIZED MASSES WILL BE COMPUTED 
INPUT GENERALIZED MASSES INSTEAD OF LUMPED MASSES 
NUMBER OF LUMPED MASSES 

TYPE OF FIT USED IN INTERPOLATING THE GENERALIZED 
AERODYNAMIC FORCE, Q(k) FOR THE INDEPENDENT VARIABLE K. 
KFIT =1 LINEAR 
KFIT =2 QUADRATIC 
KFIT =3 CUBIC SPLINE 

USE S-PLANE CURVE FIT TO AERODYNAMIC FORCES 
ASSUMING COEFFICIENTS HAVE BEEN PUT ON TAPE9 
PREVIOUSLY 

COMPUTE COEFFICIENTS FOR FIT BY CALLING SPLANE 
FUNCTION USED FOR FIT IS 

SUM(I=1,NP0LYC) OF A(I)*S**(I-1) 

PLUS 

SUM(I=1,N) OF A(NPOLYC+N)*S/(BN(l)+S) 

NUMBER OF COEFFS DESIRED IN S-PLANE FUNCTION 
= MAX(NPOLYC+N) 

IS THE CONSTANT IN THE TERM S/(BN(I)+S) 

IS THE DEGREE OF POLYNOMIAL IN S-PLANE FIT 
FOR J-TH(PRESSURE)MODE 
(DEFAULT) ALLOWS N-TH COEF FOR J-TH MODE 
TO BE COMPUTED 

(INPUT) FORCES N-TH COEF OF J-TH MODE TO BE 0.0 

Definition of Control System Model 

COMPUTE SENSOR DEFLECTION USING OVERLAY(IO.O) 

NUMBER OF SENSORS INCLUDING 2*NUMBER OF GUST VANES 
FOR ALL GUST VANES, THE FIRST NGV ARE LINEAR DISP. 

THE SECOND NGV ARE ANGULAR DISPLACEMENT 
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INTERP = 1 
ISYM =1 
NGV 

ELVN(I) 

IVANE(I) =1 
=0 

ICSREAD =1 


COMPUTE SPLINE COEFFS FOR 0VERLAY{10.0)-SENS0RS 

USING 0VERLAY(7.0) - INTERP 

MODE SHAPE SYMMETRY CODE FOR SPLINE COEFFICIENTS 

(IF 0, UNSYMMETRIC, IF 1, SYMMETRIC, IF -1, ASYMMETRIC) 

NUMBER OF GUST VANES 

DISTANCE FROM GUST VANE TO GUST REFERENCE POINT FOR 
I-TH GUST INPUT 

THERE IS A GUST VANE FOR I-TH GUST INPUT 
THERE IS NO CORRESPONDING GUST VANE 
READ IN ALTERNATE SENSOR DEFLECTION COEFFICIENTS 
TO REPLACE SOME OF THOSE COMPUTED BY SENSOR OVERLAY 
AND SAVED ON DATA COMPLEX 


Options for modifying mathematical models 


ISELECT =1 


ICHANG1=1 
ICHANG2=1 
ICHANG3=1 
ICSACT =1 

=0 

IFORCE =1 
NFORCE(I) 


IF DIFFERENT MODAL COMBINATION IS TO BE SELECTED 
FOR ANALYSIS IN CURRENT CASE THAN IN PREVIOUS NOTE: 

IF=1, THEN ALL INPUT (INCLUDING CONTROL SYSTEM) IS BASED 
ON ORIGIANAL VALUES OF NM,NC, 

AND NR. 

IF=0 AND KASE .GT. 1,NEW NM,NC,AND NR VALUE MUST 
NOT BE INPUT IN NAMELIST /INPUT/AND ALL 
CONTROL SYSTEM CHANGES SHOULD BE WITH 
RESPECT TO NM,NC. AND NR OR PREVIOUS CASE 
IF IC0NSYS=1 AND CHANGE IS TO BE MADE IN CONTROL 
NAMELIST/CONSYM/ 

IF IC0NSYS=1 AND CHANGE IS TO BE MADE IN CONTROL 
NAMELIST/FILTIN/ 

IF IC0NSYS=1 AND CHANGE IS TO BE MADE IN CONTROL 
NAMELIST/ACTINP/ 

IF ACTUATOR AND CONTROL SURFACE DYNAMIC CHARACTER- 
ISTICS ARE COMBINED AND HINGE MOMENTS AND 
INERTIAL MOMENTS OTHER THAN THE CONTROL SURFACE 
ARE NEGLECTED 

IF EFFECT OF INERTIAL AND AERO. HINGE MOMENTS UPON 
CONTROL SURFACE ARE CONSIDERED EXPLICITLY 
REPLACE SOME COLUMN IN GUST INPUT ARRAY WITH A 
COLUMN FROM AERO ARRAY 

COLUMN NUMBER IN AERO ARRAY TO REPLACE I-TH 
COLUMN OF GUST INPUT ARRAY 


Stability Analysis, PROOT 


IFLUT =1 
KVAR 

=1 

=2 

=3 

NV 


PERFORM FLUTTER ANALYSIS 

TELLS WHICH TYPE OF VARIATION IS TO BE USED 

VELOCITY VARIATION 

ALTITUDE VARIATION, DENSITY AND VELOCITY BASED ON STANDARD 

ATMOSPHERE TABLES 

DENSITY VARIATION 

NUMBER OF STEPWISE VARIATIONS 
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DV 

VO 

HO 

RHOO 

NIT 


AMOUNT OF STEPWISE CHANGE 

INITIAL VELOCITY, FOR DVAR=1 AND KVAR=3 

INITIAL ALTITUDE, FOR KVAR=2 

INITIAL DESNITY FOR KVAR=3 

MAXIMUM NUMBER OF ITERATIONS FOR SOLVING DET=0 


NCUT 

NFINE 

EPSI 

lOPTl =1 


IOPT2 =1 
ITRACE(I)=1 
=0 

PII(J) 

PII(J) 

IPS 

=0 


=1 

DEN 

lOPT =1 

XMACH 

CONFACl 

C0NFAC2 

C0NFAC3 

IGAIN =1 

DELGAIN(I,J) 


NUMBER OF TIMES STEP SIZE IS HALVED IN TRYING TO 
CONVERGE ON ROOT 

NUMBER OF DIVISIONS OF DV FOR FINE SCAN 
CONVERGENCE CRITERION 

COMPUTES INITIAL GUESSES USING 2*PIE*C/V0*0MEGA 
FOR ELASTIC MODES ONLY 

FS AND GS ARE SCALING VALUES USED TO GET REAL AND 
IMAG PARTS OF PI AND P2 

COMPUTES INITIAL INITIAL GUESSES USING STABCAR 
(DEFAULT) TRACE I-THE MODE IN PKFLUT 
DO NOT TRACE 

ARE TWO INITIAL GUESSES FOR ROOT FOR J-TH 
MODE 

IS A PRAMETER WHICH INDICATES WHICH PlI AND P2I 
ARRAYS TO USE FOR INITIAL GUESSES FOR CURRENT RUN 
P1SAVE(I)=P1I(I) AND P2SAVE(I)=P2I(I) ARE SAVED 
WHERE PI I AND P2I ARE THE FINAL VALUES OF PI AND 
P2 IN PKFLUT FROM PRECEEDING RUN OR NEW VALUES 
BEING INPUT IN NAMELIST /INPUT/ 

P1I=P1SAVE AND P2I=P2SAVE, WHICH ARE THE INITIAL 
VALUES OF PI AND P2 ARRAYS FROM PRECEEDING RUN 
INPUT DENSITY FOR KVAR=1 
MAKES DENSITY CORRESPOND TO REAL ATMOSPHERE 
ALSO INPUT XMACH, HO, CONFACl, C0NFAC2,C0NFAC3 
MACH NUMBER, USED TO COMPUTE VELOCITY FOR KVAR=2 
CONVERSTION FACTOR TO CONVERT MAIN PROGRAM VEL 
TO FT/ SEC 

CONVERSTION FACTOR TO CONVERT HO TO FT 
CONVERSTION FACTOR TO CONVERT DENSITY(IN 
SLUGS/ (FT**3))T0 UNITS DESIRED IN MAIN PROG. 
PERFORM A GAIN VARIATION 
MAKE DV=0 FOR FLUTTER OVERLAY 
AMOUNT OF CONSTANT GAIN INCREASE FOR EACH 
ITERATION FOR I-TH CONTROL AND J-TH 
SENSOR COMBINATION 


Stability Analysis STABCAR 


ISTABCR =1 
=2 

IZEROES =1 
SS(I) 


COMPUTE STABILITY CHARACTERISTICS-POLES 
COMPUTE BOTH STABILITY CHAR. AND GUST RESPONSE 
COMPUTE ZEROS ALSO 

IS AN ARRAY OF COMPLEX GUESSES FOR THE POLES 
IN STABCAR IF THE ITERATION SCHEME DOESNT CONVERGE 
ON KNOWN ROOTS. IF SS(I) IS (0.,0.),THEN THE INITIAL 
GUESS WILL BE BASED ON NEXT HIGHEST EIGENVALUE FROM 
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OMR 

THE ONE JUST CONVERGED ON 

IS AN INITIAL STARTING POINT FOR RKl AND OMEGA 

IGUST =1 

THIS MUST BE INPUT. 
RK1=0MR*C/(2.*U) 

OMEGA=CMPLX (- . 5*0MR,0MR) 

Response To Random Input (gust) 

PERFORM GUST RESPONSE ANALYSIS 

NG 

NUMBER OF GUST INPUTS USED 

NGR 

NUMBER OF GUST INPUTS READ IN 

RHO 

FOR PSD OUTPUT 
DENSITY FOR GUST NALYSIS 

NINT 

NUMBER OF INTERVALS/10 USED FOR INTEGRATION 

RFLOW 

LOWER LIMIT FOR INTEGRATION (REDUCED FREQ) 

RFCUT 

UPPER LIMIT FOR INTEGRATION (REDUCED FREQ) 

GL 

CHARACTERISTIC GUST LENGTH 

UU 

ERATION LOCATION 
VELOCITIES FOR GUST ANALYSIS 

NVEL 

NUMBER OF VELOCITIES 

lUSLF =1 

USE AN UNSTEADY LIFT FUNCTION APPROXIMATION= 

IS0UT(I)=1 

l./(l.+2.*PIE*K) 

OBTAIN OUTPUT RESPONSE TO SENSOR I 

NLOAD 

NUMBER OF LOADS 


Response To Discrete Input 


ITHRF =1 
MINT 

TMAX 

NOTE: 


IPR =1 
KFUN(I)= 


PERFORM TIME HISTORY RESPONSE ANALYSIS 

POWER OF 2 USED TO DETERMINE NUMBER OF POINTS IN 

FOURIER TRANSFORM, NPTS=2**MINT 

UPPER LIMIT ON TIME INTERVAL FOR TIME HISTORY 

MINT AND TMAX SHOULD BE CHOSEN AS FOLLOWS 

1. TMAX SUFFICIENTLY LONG TO PERMIT ALMOST STEADY 
STATE RESPONSE 

2. NPTS CHOSEN SUCH THAT MAXIMUM OMEGA IS GREATER 
THAN THE HIGHEST SYSTEM FREQUENCY OF CONCERN 

3. DELTA-T=TMAX/NPTS MUST BE ADEQUATELY SMALL TO 
DESCRIBE THE FORCING FUNCTION 

INTERMEDIATE PRINTOUR IN TIME OVERLAY 

THE NUMBER CORRESPONDING TO THE BUILT-IN FORCING 

FUNCTION DESIRED, FOR 1=1 TO NFOR 

INPUT A(I),B1(I),B2(I),B3(I),NGTH(I) - DEFINED AS 

FOLLOWS: (FOR FREQUENCY FUNCTIONS, THE TRANSFORM 

FOLLOWS; 

TYPE 1: F{T)=A*(l.-COS(Bl*T)) ONE PERIOD 

TYPE 2: F(T)=A*SIN(B1*T) 

TYPE 3: STEP IN AT 0. ,STEP OUT AT B1 

F(T)=A ,0 < T < B1 
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=0 ,B1<T 

TYPE 4: RAMP-IN , HOLD , RAMP-OUT 

F(T)=A*T/B1, 0< T < B1 

=A , B1 < T < B2 

=A* (1 . - (T-B2 )/ (B3-B2 ) ) , B2<T<B3 
-0 , B3 < T 

=0 INDICATES TIME HISTORY ARRAY INPUT TO DESCRIBE 
FORCING FUNCTION 

NUMBER OF POINTS IN TIME AND TIME RESPONSE ARRAYS 
FOR FREQUENCY RESPONSE OUTPUT TO MULTIPLE INPUT 

1 DO ONLY FREQUENCY RESPONSE 

2 DO BOTH FREQ. RESPONSE AND TIME HISTORY 

0 DO NOT DO FREQUENCY RESPONSE (ONLY TIME HISTORY) 

(SEE ALSO KFUN) 

NUMBER OF PILOT INPUTS 


Instructions For Obtaining Printed Or Plotted Results 


INTERMEDIATE PRINTOUT IN GUST OVERLAY 

SCALE FACTOR FOR TEKTRONIX PLOT. DEFAULT= .014 

INCREASE VALUE TO MAKE PLOT SMALLER 

PLOT AERODYNAMIC FORCES VS. REDUCED FREQUENCY 

PLOT FREQUENCY RESPONSE FROM GUST ANALYSIS 

PLOT POWER SPECTRA FROM GUST ANALYSIS 

PLOT TWO*GAMMA AND FREQ FROM FLUTTER ANALYSIS VS. 

VELOCITY, DENSITY, OR ALTITUDE VARIATION 

(AND DYANMIC PRESSURE IF IMATCH = 1) 

PLOT ROOT-LOCUS WITH RESPECT TO VELOCITY, DENSITY 
OR ALTITUDE VARIATION 
PLOT TIME HISTORY RESPONSE 


TABLE 4c . 

Supplementary Namelists 

DATA 10 EXPLANATION OF PARAMETERS IN NAMELIST SELECT 

IF ISELECT IS ON (ISELECT=1) 

NMODES Total number of modes to be selected from original data 

NOC(I) Mode number to be used in place of Mode I 

Example 

N0C(1)=3,N0C(2)=7,N0C(3)=9 and NMODES = 3 
Means 3 MODES will selected and they will be 
MODES 3,7, and 9 

NFNEW New NF (Number of flutter modes looked at) 

NCNEW New NC (Number of control modes) 

NRNEW New NR (Number of rigid body modes) 

NOT NECESSARY, BUT USE BOTH OR NEITHER 


IPRINT =1 
TFAC 

lAPLT =1 
IFRFPLT=1 
ISPLT =1 
IPKPLT =1 


=-l 

ITHPLT =1 


NGTH 

IFRQOUT = 
NCP 
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NKNEW 

NOK(I) 


New NK (Number of reduced frequencies) 

Frequency number to be used in place of I-TH frequency 


Sensor Location 

DATA 5 

IF ISENSE IS TURNED ON NAMELIST SENLOC MUST BE SUPPLIED 

NAMELIST /SENLOC/ XS,YS,ITYPE 
DIMENSION XS(10),YS(10),ITYPE(10) 

EXPLANATION OF PARAMETERS IN NAMELIST /SENLOC/ 

XSjYS Array of sensor locations 

for beam type, XS and YS should be unrotated 
ITYPE Indicates whether linear or angular measurements 

are to be taken 
for I = 1, NS + NS 
ITYPE(I) =1 Linear measure 

ITYPE(I) =2 Angular measure 

NSS(I) Structural section number which corresponds to I-TH 

sensor 

IF IC0NSYS=1 SUPPLY NAMELISTS CONSYM, ACTINP, AND FILTIN FOR THE FIRST CASE, 
THEREAFTER, SUPPLY CHANGES ONLY IF ICHANGl, ICHANG3 OR ICHANG2 ARE TURNED ON. 
DATA 7 

NAMELIST /CONSYM/ ASD,XKS,IORD-ISDYN, GN 

DIMENSION ASD(3,10),XKS(10),I0RD(10),ISDYN(10),GN(10(10) 

CONSYM is the block of parameters and coeffs used in subroutine SENDYN 

EXPLANATION OF PARAMETERS IN NAMELIST CONSYM 

ASD(I,J) I-TH coeff. of denominator polynomial of sensor 

transfer function which models sensor dynamics 
for J-TH sensor. 

(Default: ASD(1,J)-1. ,ASD(2,J)=ASD(3,J)=0. ) 

POLY = SUM(I=1,3) of ASD(I,J)*S**(I-1) 

XKS(J) Numerator coefficient of J-TH sensor 

(Default. XKS=1.) 

POLY = XKS(J)*S**(IORD(J)) 

lORD(J) Is the J-TH sensor type 

=0 Displacement sensor 

=1 Rate sensor 
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=2 


Acceleration sensor 


ISDYN(J) 


= 1 
=0 


GN(I,J) 


(Default = 0) 

Indicates sensor dynamics 

include sensor dynamics 

perfect sensor 
(Default. = 0) 

Gain for I-TH control and J-TH sensor 
(Default: = 0. ) 


Feedback Filter Description 


DATA 8 

NAMELIST /FILTIN/ TAU2 ,WN1 ,SN1 ,WN2,SN2 ,KIBD0,KIBD1 ,KIBD2,AN1 ,AD1 ,AL22,AL21 , 
IFILTER,AFN,AFD 

DIMENSION IFILTER(12),TAU2(12),WN1(12),SNI(12),WN2(12),SN2(12) 

C, KIBD0(12),KIBD1(12),KIBD2(12),AN1(12),AD1(12),AL22(12),AL21(12) 

C, AFN(6,12),AFD(6,12) 

REAL KIBD0,KIDB1,KIBD2 

FILTIN is the block of coefficients used in subroutine FILDYN to model filter 
logic, there is a new set of coefficients for each sensor and control 
combination. Default values for all coeffs. except WN1,WN2, and AFD(l) are 
0.. WNl and WN2 are l.OE+20 and AFD(1)=1. 


WN1,SN1 

TAU2 

WN2,SN2 

KIBDO 

KIBDl 

KIBD2 

ANl.ADl 

AL22,AL21 

AFN.AFD 


IFILTER 


EXPLANATION OF PARAMETERS IN NAMELIST /FILTIN/ 

(THESE COEFFICIENTS ARE ALL SUBSCRIPTED PER FILTER) 

Natural frequency and damping ratio at first break 
point (beginning of notch) 

Time constant at second break point (low point of notch) 
Natural frequency and damping ratio at third break point 
(end of notch) 

Gains for integral 
Proportional 

Derivative outputs, respectively 
Frequencies of zero and pole in lead-lag (lag-lead) 
compensator 

Denominator coeffs, in second order lag filter 
l./(l. + AL21*S + AL22*S**2) 

Coeffs, of numerator and denominator polynomials for 
an overall transfer function representation of control 
logic between a given sensor - control pair. 

An array indicating the types of filters to be combined 
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(in senes) 
IFILTER(I) = 


to make up a given filter 

1 No further filtering 

2 Notch filter 

3 Integral 

4 Proportional plus derivative 

5 Lead-Lag or Lag-Lead 

6 Second order lag 

7 General rational function 


Default Value is 1 

EXAMPLE: IFILTER(l )=3 ,4 , 1 MEANS FILTER WILL BE 

MADE UP OF FILTERS OF TYPES 3 and 4 


g Actuator Description 

NAMELIST /ACTINP/ AACTN,AACID, lACT 
DIMENSION AACTN(4,10),AACTD(5,10),IACT(10) 

This^array of input is for the coeffs, in subroutine ACTDYN, there should 
be NC sets of coefficients 


EXPLANATION OF PARAMETERS IN NAMELIST ACTINP 


NACT(I,J) = 
AACTN(I,J) 


0, No 

1, Yes 


AACTD 


Ith pilot command is sent to ith actuator, default off 

I-TH coefficients of numerator polynomial of actuator 
transfer function which modesl actuator for 
J-TH control . 

(Default: AACTN(1,J)=1. , All Others =0.) 

POLY = S^M(I=1,4) of AACTN(I,J)*S**(I-1) 

I-TH coefficients of denominator polynomial of transfer 
function which models actuator for J-TH control 


lACT(J) 


(Default: AACTD(1 ,J)=1. , All Others =0.) 

POLY = s2m(I = 1,5) of AACTD(I,J)*S**(I-1) 
Indicates actuator dynamics 
=1 Include actuator dyamics for J-TH control 

=0 Perfect actuator 


OUTPUT PLOTS 


NAMEL IST/PLTSEL/MODEA .MODES ,M0DEFR ,M0DEPK,M0DETH 
IF lAPLT = 1 INPUT; 

M0DEA(I,J)=1 (I,J) - TH mode curve of full aero 
Matrix including gust columns, as 
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NOTE 


NM+1 , . . . . ,NM+NG. NOTE: I and J refer to original 

values of NM. 

The NM and NC values for the parameters that follow 
are those given by select option if used. 

IF IFRFPLT=1 input; 

M0DEFR(I,J)=1 Plot (J-l)ST Derivative --J=l, Displacement 

J=2, Velocity, J=3, Acceleration--of the I"TH 
Mode Frequency Response. 

IF ISPLT=1 input; 

MODES(I) =1 Plot I-TH POWER SPECTRA. I computed as follows: 

1= 1 - INPUT POWER SPECTRA 

2, NLOAD+1 - LOADS 

NLOAD+3, NSO+NLOAD+1 - SENSORS 

NS0+NL0AD+2,NS0+NL0AD+2NC+1 - CONTROLS 

IF IPKLPLT=1 input; 

M0DEPK(I)=I Plot I-TH Flutter Mode Output 
IF ITHPLT=1 input; 

M0DETH(I)=1 Plot I-TH Time History. I computed as follows: 

I=1,NM-NC - QOFK(I) 

I=NM-NC+1,NM - CONTROL DISPLACEMENTS 

I=NM+1,NM - CONTROL RATES 

I=NM+NC+1,NM+NC+NS0 - SENSORS 

I=NM+NC+NS0+1,NM+NC+NS0+NL0AD - LOADS 

NOTE: NSO IS EQUAL TO THE SUM OF THE ONES IN THE NSOUT ARRAY 


Table 4d. 


Dynares - Order Of Input 
File: Tape 5 


Data 1 (if DCM AERO-1) 
Data 2 (if DCM AERO-1) 

Data 3 (if DCM AERO-1) 

Data 4 (if IMAT = 1) 
Data 5 (if IMAT = 1) 
Data 6 (if IMAT = 1) 


Description 


Format 


(RFRQ(I),I=1,NK) 8F10.0 

For K = 1,NK 

((AR(K,I,J),AI(K,I,J),I=1,NM),J=1,NM) 2E15.6 

For K = 1,NK 

((AGR(K,I,J),AGI(K,I,0),I=1,NM),0=1,NGR)2E15.6 
(WMASS(iLi = 1,NMAS 8F10.0 

(INERTIA(I),I=1,NMASS) 5E15.7 

For J = 1,NM 

(M0DE(I,J),I=1,NMASS) 8F10.0 

(SL0PE(I,J),I=1,NMASS) 8F10.0 

Data 4 - Data 6 Replaced by the following if: 

IMAT = 0 and DCM GMASS -1. 
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For I = 1,NM 

(GMASS(I,J),J=1,NM) 5E15.7 

Data 7 (if DCM FREQ-1) (FREQ(I),I=1,NM) 5E15.7 

Data 8 (if DCM DAMPINGS-1) (DAMPINGS(I),I=1,NM) 5E15.7 

Data 9 (if NLOADS/O and ( (CL0AD(J, I) ,J=1 ,NM) , 1=1 ,NL0AD) FREE FIELD 

DCM LOADS- 1) 

Data 10 For each IG=1,NF0R such 
that KFUN(IG)=0 

(GTH(I),TIME(I),I=1,NGTH(IG)) 2E15.6 

Data 11 (If DCM SENDEF-1) (CS(J , I ) ,J=1 ,NMNC) , 1=1 ,NS) 8F10.0 


Table 4e. 


DYNARES - TAPE 5 - DEFINITIONS 


RFRQ 

AR(K,I,J) 

AI(K,I,J) 

AGR(K,I,J) 


WMASS(I) 

INERTIA(I) 

MODE(I) 

SLOPE(I) 


GMASS(I,J) 

FREQ(J) 


DAMPING(J) 

CL0AD(I,J) 


GTH(I),TIME(I) 
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Reduced frequencies associated with Generalized 
Aerodynamic Forces. 

Real and Imaginary parts of Generalized Aerodynamic Forces 
For I-th MODE deflection and J-th mode downwash. The 
aerodynamic forces due to airplane motion are input column- 
wise. 

Real and Imaginary parts of external aerodynamic forces 
which may describe either gust force or any other external 
force such as a control surface force resulting from pilot 
commands which are not part of the automatic feedback 
system. 

Lumped masses and inertia of I-th node. 


Deflection and slope at I-th node for the J-th assumed mode. 
The mode shape data can be taken from prior vibration 
analyses or represented by some suitable mathematical 
function. 

Generalized mass . 

Modal frequency from the vibration analysis or term used 
to describe the modal stiffness. 

Modal structural damping ratios. 

Stress coefficients at the I-th station for the J-th mode 
used in equation . The terms though meaningful only 
for slender beams, can be used to define any stress in 
any type of structure for which modal analysis is 
appl i cable. 

An array which describes an external forcing function. 

GTH(I) is the amplitude at some time TIME(I). Used in 
the time history overlay. 



DYNARES - Sample Input 
Control Deck 


[>YNA3»T5O0. BINXX XXXXX 

USERi XXXXXXXt XXl.'XXXX. 

CMARGEiXXXXXXiLRC. 

HAP, OFF. 

GET, IHITIAL/UH-887010C. 

CALL,INITTAL. 

GET, DYNAU18/UN»G87010C. 

GET,TAPEj^=BIS32. 

GET,Tf.?E2=3ISC. 

ATTACH, FTNHLIB,LRCGOSF/UH«*LIBRARY. 

GET,ACPOLIB/UN=887018C. 

LDSET,LIB=ACPOLIB/FTNHLIB/LRCGOSF. 

L0AD<DYHAU1B> 

DYNARES. 

PL0T.UARIAH<FSH»8. ,FSU*11. I AUTO(0)> 

END OF FILE . 


In the above run it is presupposed that all of configuration data has been 
stored on the file B1S92. If this is not the case, these data may be input 
through TAPES compiled as instructed in table 4c. 
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ISAC - Sample Problem for consecutive run 


ISACB,T0550. 

USER, 

CHARGE, 101390, LRC. 

MAP, OFF. 

GcT,IHITIAL/UH»887010C. 

CALL, INITIAL. 

GET,TAPE5=BIS5. 

GET,DCf1/UH=887010C. 

DCH. 

REPLACE, TAPE9=BIS92. 

RETURN, TAPES. 

GET,TAPE1=NASBISC. 

GET,TAPE3=DLIBC. 

GET,TAPE5=DLhBC. 

GET , LGO=DL IHBN/UN=8870 1 0C . 

LDSET<LIB=ACP0LIB/FTNMLIB^LRCG0SF) 

LGO. 

REPLACE, TAPE9=BIS92. 

DAYFILE, ISACDAY. 

REPLACE, ISACDAY. 

PLOT.UARIAN<FSH=ll. ,FSU=8. ,AUTO<0)) 
REWIND, TAPES. 

RETURN, TAPEl. 

GET,LGO=DLATBH/UN=887010C. 

LDSET<LIB=ACPOLIB/FTN11LIB/LRCGOSF) 

LGO. 

REPLACE, TAPE9=BIS92. 

DAYFILE, ISACDAY. 

REPLACE, ISACDAY. 


STATUS, F. 

REWIND, OUTPUT. 

COPYEI,OUTPUT=AA. 

REPLACE, AA. 

RETURN, TAPEl. 

ATTACH, FTNMLIB.LRCGOSF/UN-'tLIBRARY. 

GET,DYLIB,ACPOLIB/'UN*837010C. 

GET,TAPE9=BIS92. 

GET,TAPE2^BISC. 

GET , LGO=DYNAUlB/UNs88701OC. 
LDSET<LIB=ACP0LIB/FTNMHB/LRCG0SF> 
LGO . 

PL0T.UARIAN<FSH=8. ,FSU=11. ,AUTO<0)) 
DAYFILE, ISACDAY. 

REPLACE, ISACDAY. 

EXIT. 

DAYFILE, ISACDAY. 

REPLACE, ISACDAY. 

STORE 

GHASS -1141 
FREQ -1121 
DAMPINGS -1121 
END 

(2E15.7) 

BISP C 

SAME 

BISP C 

SAME 

BISP C 

TOC 

END 


NOTES 


1. TAPEl is a formatted tape generated by the vibration analysis 

program. 

2. On an initial run. 

3. Input data and results are shown in table 5 and figures A2 

through A4. 
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FIGURE 1. BLOCK DIAGRAM BASIC AIRPLANE, AUTOMATIC CONTROLS 
PILOT INPUTS AND GUST INPUTS 
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Figure 2. Program Structure 















Blsplinghoff jet trans-cant 


FIGURE A2. STRUCTURAL NODE AND AERODYNAMIC BOXES 



Bisplinghoff jet trans-cant Mode ' 1, Freq - 1. 83HZ, Gmass ■ 897E+01 



Bisplinghoff jet trans-cant Mode - 2 , Freq - 3. 65HZ, Gmass - . 995+03 


FIGURE A3. VIBRATION MODE SHAPES 
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CTITU • i ISP- JET TRAMSPORT 
«$U»TXTLC* CANTUEWER 


tlACEt • 
CEICCMUECTOR 
CRCAl OUTPUT 
•SUICASE ID • 


1 

CEICENUALUE • 10 

.1322eilE*B3 RODE • 

1 

G 

a.E^eo 

-COHT- 


B.E4O0 

6 

G 

0.E>«O 

-CONT- 


8.422469E-04 

12 

C 

e.ET«« 

-CONT- 


l.058fi48E-#3 

17 

G 

B.E*M 

-COKT- 


2.673194E-83 

22 

G 

B.E^M 

-CORT- 


3.22»«44E-a3 

27 

G 


-CONT- 


3.377827E-a3 

fTITLE - BISP- 

JET 

TRANSPORT 


•SU>TITl£« CA«TILEVER 

«LP»EL 

f£IC£Jft«CTO« 

•REAL OUTPUT 

tSUtCASE ID • 2 

SEICEMUftLUE ■ 60.E27E27SE+R3 «0Dt • 


-CONT- 

1 

C 

O.E«^B3 

B.E-fBG 

-CONT- 

e 

G 

6 E+e« 

-S.584B62C-*': 

-CONT- 

12 

C 

B.e^ 

I.SIOSB-IC-B^ 

-OONT- 

17 

G 

e.E^OB 

2.ieiaeec-S3 

-CONT- 

22 

G 

• .E>M 

5.S3SS83E-C3 

-CONT- 

27 

C 

B.E-»«3 

6.4H23EC-B3 

rriTLE 

• BISP- 

-JET TRANSPORT 

aSUBTITU* CAMTIUUER 
aUUEL • 
•EIGE31UECTOR 
aSEAt. OUTPUT 


•SUBCASE ID • 
7 


3 


1 

3 

6 

7 


a 

8 




10 

B.EtBB 

e.E«ee 

11 

f.t*99 

e.E*^e9 

12 

9.1*99 

5.e4339EE-B2 

13 

3.O41B8BE-05 

a.E4^ea 

14 

«.£*«» 

1.995S29E-ei 

15 

S.931993E-B5 

9.E+B9 

16 

o.E^oe 

3.967C01E-01 

17 

4.S2S544E-B5 

0.E+00 

18 

e.E«^e« 

7.ee40isE-ei 

19 

3.482329E-0S 

a.E^oa 

20 

o.E^-eo 

i.se9aeeE+00 

21 

2.e24S96£-eS 

1 

e.E^ee 

22 

23 

25 

23 

29 

30 

31 

32 

B.E+B9 

O.E-fBO 

33 

e.E>eo 

O.E+08 

34 

• .E*M 

-3.1Hei3E-02 

35 

3.C1797CE-92 

O.C»80 

36 

B.E+09 

-1.728147E-01 

37 

C.94548SE «2 

C.E>08 

3a 

e.Ei^eo 

-9.132470E-03 

39 

S.254SeiE ca 

C.ZH9 

4o 

e.E^ee 

4.5ieoi3E-ei 

41 

•G.4452eSE-02 

0.E«83 

42 

o.E^ee 

1.000300E«0O 

43 

6.S43473E-02 

O.E*00 

44 

45 
47 

50 

51 

52 

53 


TABLE 5a. sample PROBLEM INPUT D.ATA, NASBISC 






YOR « a3E*01, 

BISPLINGHOFF JET TRANS-CAHT 

IPLAN » 1, 

NASTRAN 



IBOTH « 1, 

$DLINPT 

NSECTNS«l,NSECTHA*l,NDF«6i 

ISTOP =0, 

NNODES* 

6f 

NH0DET=6,NPMX*6, 

$END 

IPLTNOD 


,IPRINT*1, 

$MODEPLT ISTOP-1 , MODEMS-0, 1 , SC2DEF— 

NSKIPL 

,*? 

,SCX=108.|SCY=100., 

la$ 

IPLTMS= 

1, 

IBAUD=120» 

DISP -1 1 

«ODE=l, 

U 


MODES 1 1 

IDCH=1, 

IPLTBOX=i$ 

HHD 0 1 

$MODEPLT 


SPLINE 1 1 

INODELh 

1=1 

f 

END 

NSURF 

ss 

1 f 

1,8,0,278.75,0. ,0a 

XLE 

= 

.2Ei03» .244E-f03, 

1»78,1 

YLE 

s 

0a 0y a5E‘^03f 

0a, 0a, 1,3, ‘^5 

NLEP 

s 

2, 

90a ,0a ,2,3, *^5 

XTE 

s 

.425E+03, .3438E+03i 

186a ,0. ,3,3, +5 

VTE 

s 

0t0y a5E+03y 

268a ,0a ,4,3, ^5 

NTRP 

s 

2, 

368a ,0a ,5, 3, *^5 

NSPH 

s 

19, 

458.,0.,6,3,4>5 

HCHD 

s 

3, 

7 

• 

HSSP 

s 

1 1 


SCXMS 

ss 

.5E+02, 


SCYMS 

s 

•75E+02, 


SCZMS 

= 

alE+01, 


SC2DEF 

s 

.8E+00, 

, 

XU 

s 

.llE+04, 


YU 

= 

a 5E+02, 


2U 

sz 

a5Er01, 


ALPHA 

s 

a45E+02, 


BETA 

3 

-a5E+01, 


XTRANS 

B 

0a 0, 


HE 

s 

2, 


MODEMS 

=1 

,0,0,0,01 


XOR 

s 

a2E401, 




TABLE 5b. SAMPLE PROBLEM INPUT DATA, DLIBC 



BISPLINGHOFF JET TRANSPORT CANTILEUER 

♦DLATINP FMACH«. 7, ACAP-81250. , NDELT«1 , NP»1 , NB«0, NRF«6, RFREQ*0. , . 2» . 4, 

. 6, . 8, 1 . , NSTRIPalS, NSU«0, HBU«0, NM0DESs2, REFSPN^l . , NBOX-78, HSNC»98, 
REFCH0S162. 5, IPRINT-l , XZERO-200. , XCG-280. 4, NR*0, NC»0, NELAST*2, IDCM*1$ 

AERO 1 1 
HHD 1 1 
END 

200. 425. 243.8 343.8 0. 500. 

0. 0. 14 7 1. 

0. .167 .333 .5 .667 .833 1. 

0. .08 .16 .24 .32 .40 .48 .56 .64 .72 .80 .88 .96 1. 

16 7 12 13 18 19 24 25 30 31 36 37 42 43 48 49 54 55 68 61 66 67 72 
73 78 
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TABLE 5c. SAMPLE PROBLEM INPUT DATA, DLABC 



d) INPUT DATA, BISC 


1 

BISPLINGHOFF JET TRAN FLUT CANTILEUER 

♦ INPUT NM=2, NR*0, NC=0, NK=6, NLOAD=0, C*162. S, KFIT«s2, NS«0, XCG-280. 4, 
NNODES=Si NSECTNS*! , NPMX=6, IPLATE«0, 1 , NINTs30, GL«30000. , 

XMACH«. 7, H0«50000. , CONFACl*. 083333, C0NFAC2«l . , C0NFAC3=, 000048225, 
DU*300. ,NU=50, IOPT1«0, !GAIN=0, 

IFLUT=1,KUAR=1,IPKPLT«1, 

RHO0=1.15-7,IOPT2=1, 

DEN=1. 146-7, U0=10000., 

RFCUT^l . 0, NCUT=4, NFINEM, NIT=26, ISYM=1 , ICSACT*! , EPSI*! . E-6, RFLOW*0. 
ISENSE=0, IPRIHT=10,NUEL=1, ICONSYS«0,KUAR=1, ISTABCR=0, IDCH=:1, 
ISELECT=0, 

FS=.01,GS=1.,ND=1, 

R0=0. , 1 . , X0=27S. 8,0. , Y0-0. ,0. , ISS-1, 6, 8, 181 
ALL 1 0 
END 

♦PLTSEL HODEPK»l,l$ 


e) INPUT DATA, BISS 


8.9735 0.0000 

0.0000 994.5340 

1.8300 3.6540 

.0300 .0300 


TABLE 5. CONCLUDED 



End of Document 



